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Сформулирована модель формирования диффузионной зоны между материалом и покрытием в процессе импульсной элек- 
тронно-лучевой обработки. Разработан алгоритм численного решения задачи неизотермической диффузии, учитывающий раз- 
личие пространственных и временных масштабов теплопроводности и диффузии в твердых веществах. Исследована зависи- 
мость характеристик диффузионной зоны от режима облучения. 


При моделировании поверхностной термической Вследствие различия пространственных и вре- 
обработки материалов; формирования переходной менных масштабов тепловых и диффузионных 
зоны между покрытием и подложкой в процессе им- процессов в конденсированной фазе такие задачи 
пульсной электронно-лучевой или лазерной обра- оказываются весьма непростыми для их численной 
ботки систем типа «пленка-подложка» возникает реализации, так как требуют больших вычисли- 
необходимость совместного решения системы ура- тельных ресурсов или специально разработанных 
внений теплопроводности и многокомпонентной численных алгоритмов. 

диффузии, зачастую включающих источники или Имеющиеся в литературе математические мо- 
стоки тепла и массы вследствие физико-химических дели процессов многокомпонентного диффузион- 
превращений. Подобные математические модели по- ного переноса по той же причине ограничены, как 
зволяют существенно дополнить экспериментальные правило, изотермическими приближениями. Гото- 
исследования за счет подробного изучения динами- вые пакеты прикладных программ зачастую рас- 
ческой картины формирования диффузионной зоны. считаны на достаточно узкий круг задач. 
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Цель настоящей работы заключается в постро- 
ении математической модели формирования диф- 
фузионной зоны в поверхностном слое материала, 
контактирующем с покрытием, в процессе обра- 
ботки импульсным источником энергии и в ее чи- 
сленном исследовании. 


1. Математическая постановка задачи 


Предположим, что покрытие малой толщины 
/г,, предварительно нанесенное на подложку, содер- 
жит в своем составе легирующие элементы, кон- 
центрации которых обозначим как й„ /= 1,2. В на- 
чальный момент времени 1=0 легирующие элемен- 
ты распределены в покрытии равномерно, а в под- 
ложке отсутствуют. При 1>0 система со стороны по- 
крытия подвергается действию импульсно-перио- 
дического источника тепла, например потока элек- 
тронов. Полагаем, что радиус пучка электронов 
превышает поперечные размеры образца, что по- 
зволяет ограничиться одномерной постановкой за- 
дачи, которая включает в себя уравнение теплопро- 
водности и уравнения диффузии легирующих эл- 
ементов из покрытия толщины в основу (рис. 1). 

< 7(0 



Рис. 1. Иллюстрация к постановке задачи 
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- диффузионные потоки легирующих элементов, 
Д,, - парциальные диффузионные коэффициенты; 
с, р, Я - теплоемкость, плотность и теплопровод- 
ность подложки. 


Малая толщина слоя присадочного материала 
по сравнению с толщиной подложки и толщиной 
теплового пограничного слоя, формирующегося в 
подложке за время обработки, позволяет не рас- 
сматривать распределение температуры и концен- 
трации в этом слое. Тогда в качестве граничного 
условия для температуры при х=0 будет справедли- 
во граничное условие в виде дифференциального 
уравнения: 

л = 0 : Я Ц- = /?,с,р, д -^ + се (Г 4 - Д 4 )- 9 {1), (3) 

дх ді 


где Т е - температура окружающей среды, 
ае( Ѵ—Т*) - потери в окружающую среду (вакуум), 
ст - постоянная Стефана-Больцмана, е — коэффи- 
циент черноты. Вид функции ц(1) и начальные 
условия зависят от описываемой стадии технологи- 
ческого процесса. Во время действия внешнего им- 
пульсно-периодического источника имеем 

Ф) = ?о/(0; А*) = ІШ 

к = 1 

где 

IX (і,+І р )(к-])<К1 к + (( +І р )(к-1); 

X х* -!)**< (4 Н)*» 

I, - длительность импульса; і р - длительность пау- 
зы; к - номер импульса облучения; и - число им- 
пульсов. Плотность мощности, измеряемая в 
Дж/(см 2 -с), есть <?(/); энергия в импульсе - #(?)С 

Если источник не действует по истечении п им- 
пульсов, то #(/)= 0. 

Для концентраций на границе имеем условия: 

* = 0 : и ' ІЬ = и ' ІИ = “ (4) 

На бесконечном удалении от нагретой области 
источники и стоки тепла и массы отсутствуют: 

х— »оо:Я — = 0; (5) 

дх 


У, = 0; У 2 = 0. (6) 

Зависимости теплофизических свойств от тем- 
пературы аппроксимируются полиномами. Коэф- 
фициенты диффузии в общем случае зависят от 
концентраций по степенному закону, а от темпера- 
туры - по закону Аррениуса. Покрытие и подлож- 
ка имеют различные температуры плавления Т рКЬ 
Т рк . Следовательно, для описания температурных 
зависимостей свойств имеем соотношения: 

А, = (А/ + А* а 1 + А і а 2 + Оу а г а 2 + Оу а~+ АѴ) Х 

хехр(-| : ), ,7, 

Я = Я 0 + Я, Г + Я 2 Г 2 + Я3Г 3 + Я 4 Г 4 , (8) 
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Соотношение, аналогичное (9), имеет место для 
подложки. Здесь 5 есть 5 - функция Дирака. 

Начальные условия к задаче (1) - (5) имеют вид 


^ 0 а, ^10 5^2 А) * 9 

а, = 0, а 2 = 0, д > 0 


Т=Т 0 ,х> 0. (10) 


На основе модели требуется исследовать влия- 
ние параметров источника на динамику изменения 
температуры и перераспределение концентраций 
легирующих элементов. 
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2. Результаты численного исследования 

модельной задачи 

Задача решалась численно по неявной консер- 
вативной разностной схеме с использованием мето- 
да прогонки. В настоящей работе был применен 
специальный алгоритм для численного решения за- 
дач такого типа, описанный в [1]. Задачи теплопро- 
водности и диффузии решались на разных, но со- 
гласованных, разностных схемах. Шаги по про- 
странству и по времени выбирались из соотноше- 
ния характерных диффузионных и пространствен- 
ных масштабов для процессов теплопроводности и 
диффузии. В качестве характерного пространствен- 
ного масштаба для диффузии использовалась тол- 
щина диффузионного пограничного слоя, форми- 
рующегося в течение времени, равного длительно- 
сти одного импульса /. В расчете пространственно- 
го масштаба используется минимальный из коэф- 
фициентов диффузии (рассчитанный при темпера- 
туре, близкой к температуре плавления подложки): 



Пространственный масштаб для теплопроводно- 
сти может быть вычислен аналогично 



Тогда шаги разностной сетки в задачах теплопро- 
водности и диффузии связаны следующим образом 

Ах т /Ах 0 = ^/(срО тЬ ) = і/ л/Гё, 

где Ье^ицн/к: - число Льюиса для данной задачи. 
Аналогично можно выбирать и соотношение меж- 
ду шагами по времени. Но использование неявных 
абсолютно устойчивых разностных схем позволяет 
об этом тщательно не заботиться. 

Еще одна проблема при решении задач неизо- 
термического массопереноса в конденсированных 

в б 


средах заключается в следующем: известные усло- 
вия на матрицу коэффициентов диффузии, полу- 
ченные для изотермических условий и необходимые 
для обеспечения неотрицательности решения, не 
всегда выполняются для неизотермических задач и 
задач с коэффициентами диффузии, зависящими от 
температуры и концентраций. Эта проблема реша- 
ется при построении численного алгоритма исполь- 
зованием ограничений на коэффициенты модели, 
следующие из неравновесной термодинамики [2]. К 
сожалению, имеющиеся в литературе эксперимен- 
тальные данные для парциальных коэффициентов 
диффузии, полученные при проведении «изотерми- 
ческих» экспериментов, часто не удается использо- 
вать, как не удовлетворяющие этим условиям. Поэ- 
тому проблема определения парциальных диффузи- 
онных коэффициентов остается актуальной. 

В расчетах, описанных ниже в модельной задаче, 
использованы следующие параметры: 7^=1710 К, 
7^=1810 К, АН), 72-1, 37Э0- 7 Г-6, 00- ІО- 7 Р+3,39- Ю- 1 ^, 
А, =0,05 см; С Ѵ =С Ѵ |=0,68 Дж/г-К, С ь =С и = 0,65 Дж/г-К, 
р=р,і=7,61 г/см 3 , р,=рц= 7,865 г/см 3 , Г„,=250 Дж/г, 
4=400 Дж/г, ДѴ= I К)-' смУс, 7)«=Гю- 2 смУс, 
А° 2 =А°1=0, Е= 140000 Дж/г. 

Расчеты проведены для двух различных частот 

- Г = — - — = 50 и Г=10 Гц. 

{ і +1 р 

Анализ результатов численного решения задачи 
показывает, что распределение и максимальное 
значение температуры и концентрации в обраба- 
тываемом образце зависят от максимальной плот- 
ности мощности, числа и длительности импульсов, 
а так же от частоты источника. 

На рис. 2, а-в, в качестве примера представлена 
зависимость температуры от времени в точке х=0. К 
концу первого импульса (рис. 2, а) длительностью 
0,005 с температура в покрытии и на поверхности 


а 








Рис. 2. Зависимость температуры и концентрации легирующих элементов от времени в точке х=0 при следующих параметрах: 
%=0,005, І р =0,015 с, Цо =6- ІО 4 Вт/см 2 ; в, е) п=5, 6, д) п=2, а, г) п=1 
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подложки при данном наборе параметров не пре- 
вышает 1100 К. За это время концентрация леги- 
рующих элементов в покрытии практически изме- 
няется слабо (рис. 2, г). Диффузия наблюдается 
ли ть при температуре 71>700...900 К, которая су- 
ществует в поверхностном слое лишь в течении 
Д/»0,02 с. Второй импульс, действующий за первым 
после времени паузы ^,=0,00 1 5 с, приводит к боль- 
шему повышению температуры, которая ускоряет 
диффузию так, что за время второго импульса и при 
остывании до температуры I концентрация эл- 
емента в покрытии уменьшается до С/М), 0179, а 
С 2 ’*0,297. Пять импульсов (рис. 2, в, ё) достаточно 
для того, чтобы концентрация первого элемента С, 
в покрытии уменьшалась почти до нуля, а второго 
до С 2 *»0,045. Пики на кривой зависимости темпера- 
туры от времени соответствуют окончанию импуль- 
сов, плато - плавлению (во время действия третье- 
го импульса и кристаллизации - при остывании по- 
сле 5 импульсов). Концентрации легирующих эл- 
ементов в покрытии во времени практически отсле- 
живают поведение температуры (рис. 2, е). Про- 
странственное распределение температуры и кон- 
центраций в подложке в различные моменты вре- 
мени, представленные на рис. 3, иллюстрируют тот 
факт, что ширина диффузионной зоны х в , которую 
можно определить по условию 


тіп 


[ с,-С,.(/,х д ) ] 


< 0 , 01 , 


много меньше ширины зоны прогрева х т , где тем- 
пература превышала 700... 900 К в ходе обработки. 

Численное значение С', С 2 *, в покрытие к окон- 
чанию процесса обработки и величины х т , х в , суще- 
ственно зависят от параметров внешнего воздей- 
ствия. Такие зависимости представляют интерес 
для конкретных систем. 

Учет перекрестных диффузионных потоков и за- 
висимостей коэффициентов диффузии от концен- 
траций, так же оказывает существенное влияние на 
динамику процесса, что на рисунках не показано. 



Рис. 3. Пространственное ра- 
спределение температу- 
ры и легирующих эл- 
ементов в моменты вре- 
мени 0,03, 0,04, 0,05, 
0,07 при следующих па- 
раметрах: п=5, уі=0,01, 
І р =0,01 с, д 0 =3- ІО 4 Вт/см 2 

Замечено, что при уменьшении плотности мощ- 
ности источника перераспределение концентраций 
происходит более плавно, но для их существенного 
изменения требуется большее число импульсов. 



При частоте 10 Гц влияния параметров источ- 
ника на распределение концентраций легирующих 
элементов в образце (< 7 0 =Ю 4 и /,=0,01, /,,=0,09) прак- 
тически не заметны. 

3. Распределение температуры и концентраций 
в системе «медное покрытие - железная подложка» 
при импульсной электронно-лучевой обработке 

Предложенная модель была использована для 
численного исследования перераспределения кон- 
центраций в системе «медь - сталь», эксперимен- 
тальные исследования модификации которой с по- 
мощью импульсного потока электронов изложены 
в [3]. В экспериментах [3] система пленка-подлож- 
ка подвергалась импульсному плавлению пучком 
электронов с низкой энергией (20... 30 кэВ) и высо- 
ким потоком (2...3 мкс, 2... 10 Дж/см 2 ). Авторы по- 
казали, что многократное импульсное плавление 
данной системы приводит к увеличению тол щ ины 
диффузионного слоя на глубину, существенно пре- 
вышающую глубину жидкофазной диффузии, что 
значительно улучшает прочность полученного сое- 
динения. Модель ( 1) — ( 10) полностью соответству- 
ет условиям эксперимента этой работы. 

Численные исследования были проведены без учета 
наличия в стали легирующих элементов. Свойства ве- 
ществ - основы (Ре) и покрытия (Си), использованные 
в расчетах, взяты из [4, 5]: 7^=1358 К, Т ша = 1810 К, 
Я, .„,=0,45 Вт/см-К-с, /г,=0,000512см; С „„=0,68 Дж/г-К, 
С „„=0,65 Дж/г-К, р, „=7,61 г/см 3 , Рі „=7,8651 г/см 3 , 
Д*_=250 Дж/г, і рМ =20500 Дж/г, С іС =0,43 Дж/г-К, 
Сі о,=0,5 Дж/г-К, р, с „ = 0,88 г/см 3 , р ІС =0,85 г/см 3 , 
Я?, о=0,33 см 2 /с, Г=1 40000 Дж/г. 




Рис. 4. Зависимость температуры и концентрации легирующего 
элемента от времени в точке х=0 при р 0 =3,32- ІО 4 Вт/см 2 : 
а) п=2, 1;=1р=0,01 с, 6) п=4, Ь=1 Р =0,005 с 


Для данной системы характерно нелинейное 
распределение концентрации меди, что связано с 
нелинейностью задачи и учетом плавления (рис. 4). 
Результат существенно зависит от числа и длитель- 
ности импульсов облучения даже при постоянстве 
энергии, введенной в систему за время обработки. 
Так, для рис. 4, а, и рис. 4, б, плотность энергии, 
введенной в вещество, одинакова и равна 
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6,64- ІО 2 Дж/см 2 . Если эта энергия вводится за два 
импульса, концентрация в поверхностном слое (на 
границе раздела меди и железа) к моменту остыва- 
ния системы до температуры Г=900 К уменьшается 
до С*«0,74; если та же энергия вводится за 4 им- 
пульса с меньшей в два раза длительностью (так что 
полное время обработки остается тем же), то полу- 
чаем С*«0,79. Ширина диффузионной зоны, фор- 
мируемой в течение обработки и остывания систе- 
мы в первом случае составляет ^0,0059, а во вто- 
ром - _ХдЯ=0,0066 (рис. 5). Качественное распределе- 
ние концентрации легирующего элемента в покры- 
тии практически не отличается при варьировании 
параметров источника. 



Рис. 5. Пространственное распределение Си в Ре в различные 
моменты времени а) 1 ~ 0,0125; 2 - 0,025; 3 ~ 0,0375; 
4 - 0,056)1-0, 0125; 2 - 0, 025; 3~ 0, 05; 4 ~ 0,1 с параме- 
трами источника До =3,32-10' Вт/а/?: а) п=2, 1 : =1: р =0,01 с, 
х^О.0059; 6) п=4, і, =1 Р =0,005 с, х^0,0066 


Ширина зоны прогрева х т и ширина диффу- 
зионной зоны х в зависят и от плотности мощности 
д 0 (т. е. от энергии, вводимой за 1 импульс) при по- 
стоянстве вводимой энергии (рис. 6). Обработка 
одинаковым числом импульсов с различной вели- 
чиной (и соответственно с различной частотой) 
приводит к различным значениям зон прогрева и 
диффузии. Зона прогрева, определяемая по задан- 


ному значению температуры Г=900 К, увеличива- 
ется в процессе нагрева и плавления, при остыва- 
нии, начиная с некоторого момента времени, - 
уменьшается (рис. 6, б). Диффузионная зона меня- 
ется (увеличивается), пока температура достаточна 
для диффузии, а затем перестает изменяться (рис. 
6, б). Оказалось, что даже незначительное увеличе- 
ние частоты внешнего воздействия с Г =(? і +/ ; ,) _1 до 
Г= 100 Гц приводит к существенному увеличению 
итоговой ширины диффузионной зоны (соответ- 
ственно кривые 2 и 1 на рис. 6), что качественно со- 
ответствует экспериментальным данным [3]. 



а б 

Рис. 6 Ширина зоны прогрева а) и ширина диффузионной 
зоны б) в процессе обработки материала с параме- 
трами источника: 1) ц а =3, 32104 Вт/см 2 , іі=І р =0,005 с, 
п=4; 2) Цо=2,5- ІО 4 Вт/см 2 , Ь=і р =0, 00644 с, п=4 

Таким образом, в работе предложена модель 
формирования диффузионной зоны в процессе об- 
работки материала с покрытием внешним источ- 
ником энергии на основе экспериментальных дан- 
ных. Показано, что режим облучения существенно 
влияет на распределение концентрации элементов 
и тем самым на формирование переходной зоны 
между покрытием и подложкой. При расчетах для 
конкретной системы получены результаты, каче- 
ственно соответствующие наблюдаемым законо- 
мерностям. 
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